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are |ftf , 30 | < 0.23, |/if , 40 | < 0.020, |h7 0i30 | < 0.23, and |ft^ 0>40 | < 0.019 for A = 1000 GeV. 
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Studies of events containing pairs of vector bosons pro- 
vide important tests of the standard model (SM) of elec- 
troweak interactions. In the SM, the trilinear gauge cou- 
plings of the Z boson to the photon are zero; there- 



fore, photons do not interact with Z bosons at lowest 
order. Evidence for such an interaction would indicate 
new physics 0, ■ 

Studies of Z boson and photon production have been 
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made by the CDF jjj and D0 collaborations using pp 
collisions, and by the DELPHI 0, L3 @, and OPAL 
collaborations using e + e~ collisions. We present a new 
study of Z7 production using Z boson decays to e + e~ 
and At + /i~, where the dilepton system can be produced 
by either an on-shell Z boson, or a virtual Z boson or 7 
(the Drell-Yan process). The dilepton plus photon final 
state, £ + £~j, can be produced in the SM through either 
of two processes. The photon may be emitted through 
initial state radiation (ISR) from one of the partons in 
the p or p, or produced as final state radiation (FSR) 
from one of the final state leptons. We collectively refer 
to these processes as Zj production. 

The SM Zj processes produce photons with a rapidly 
falling transverse energy, E T . In contrast, anomalous 
ZZj and Z77 couplings, which appear in extensions 
of the SM, can cause production of photons with high 
E T and can increase the £ + £~^f cross section compared 
to the SM prediction. Below we describe a search 
for this anomalous production within the framework of 
Ref. ||- This formalism assumes only that the ZV7 
(V=Z, 7) couplings are Lorentz- and gauge-invariant. 
The most general ZV"f coupling is parameterized by two 
CP-violating {h\ and h%) and two CP-conserving (h^ 
and h\) complex coupling parameters. Partial wave uni- 
tarity is ensured at high energies by using a form-factor 
ansatz = hj /{l + J/A 2 )™* (i = 1,..., 4), where Vs 
is the parton center-of-mass energy, A is the form-factor 
scale, and rii is the form factor power. We set the form 
factor powers n\ = 713 = 3 and 77,2 = n± = 4, in accor- 
dance with 

The data are collected by the D0 Run II detector at 
the Fermilab Tevatron Collider with pp center-of-mass 
energy ^/s = 1.96 TeV between April 2002 and June 
2004. The integrated luminosities used for this analy- 
sis are 320 pb^ 1 for the electron final state and 290 pb _1 
for the muon final state. 

The D0 detector || consists of an inner tracker, sur- 
rounded by liquid-argon/uranium calorimeters, and a 
muon spectrometer. The detector sub-systems provide 
measurements over the full range of azimuthal angle <p 
and over different, overlapping regions of detector pseu- 
dorapidity r\. The inner tracker consists of a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both located within a 2 T superconducting solenoidal 
magnet. The CFT and the SMT have coverage out to 
\r\\ < 1.8 and |?7| « 3.0, respectively. The calorimeter 
is divided into a central calorimeter (CC) |r;| < 1.1 and 
two end calorimeters (EC) housed in separate cryostats 
which extend coverage to |?7| ~ 4. The calorimeters are 
longitudinally segmented into electromagnetic (EM) and 
hadronic sections. The muon system lies outside the 
calorimeters and consists of tracking detectors, scintil- 
lation trigger counters, and a 1.8 T toroid magnet. It 
has coverage up to |ry| « 2.0. Luminosity is measured us- 
ing plastic scintillator arrays located in front of the EC 



cryostats, covering 2.7 < 1^1 < 4.4. 

The data are collected with a three-level trigger sys- 
tem (LI, L2, and L3). We require that the events in the 
electron decay channel satisfy one of the high-i?T sin- 
gle electron triggers, while the events in the muon decay 
channel must fire one of the high-py single or dimuon 
triggers. The single electron triggers require a signifi- 
cant amount of energy deposited in the EM section of 
the calorimeter at LI. At L3, additional requirements 
are imposed on the fraction of energy deposited in the 
EM calorimeter and the shape of the energy deposition. 
The efficiency of the electron trigger requirement is about 
80% for an electron with Et ~ 25 GeV and more than 
98% for Et > 30 GeV. The muon trigger requires hits in 
the muon system scintillator at LI, and in portions of the 
data set also requires spatially-matched hits in the muon 
tracking detectors. At L2, muon track segments are re- 
constructed and pr requirements are imposed. At L3, 
some of the triggers used in this analysis require muon 
candidate events to have a high-py track reconstructed 
in the inner tracker. The logical OR of single and dimuon 
triggers has an efficiency of 92% for muons from Z boson 
decay. 

Electrons are reconstructed as clusters of energy in the 
calorimeter. These clusters are required to have 90% of 
their energy deposited in the EM calorimeter (in either 
the central calorimeter \r]\ < 1.1, or the end calorimeter 
1.5 < |?7| < 2.5). We require that the longitudinal and 
transverse shower shape of the cluster is consistent with 
that expected from an electron, and that the cluster is 
isolated from other activity in the calorimeter. Electron 
candidates in the central calorimeter are required to have 
spatially matched tracks. At least one electron candidate 
must be identified in the CC region and at least one is 
required to have pr > 25 GeV/c. Muons are identified 
by a central track matched to segments in the muon sys- 
tem. The muon must be within |?7| < 2.0. To reduce 
potential contamination from hadronic bb events, we im- 
pose isolation requirements on the muon candidates in 
both the calorimeter and central tracker. To remove the 
background from cosmic ray muons, muon tracks must 
originate from the beam region and not be back-to-back. 
Z boson candidates are reconstructed by requiring a pair 
of high-px {pt > 15 GeV/c) electrons or muons that form 
an invariant mass above 30 GeV/c 2 . 

In addition to a Z boson candidate, we require events 
to have a photon candidate, with a separation from both 
of the leptons of ATI = y/(</> e - (f) 1 ) 2 + {m - ry 7 ) 2 > 0.7 
and with EZ > 8 GeV. Photons are reconstructed as en- 
ergy clusters in the central calorimeter. The transverse 
shower shape of the cluster must be consistent with that 
expected from a photon. We also require a photon can- 
didate to deposit at least 90% of its energy in the EM 
calorimeter and to be isolated from other activity in the 
calorimeter and the tracker. 

Muon and electron detection efficiencies for the above 
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requirements are determined using a sample of Z — > ££ 
events. In the electron channel the combined trigger and 
reconstruction efficiency is measured to be (73 ± 4)%. In 
the muon channel it is measured to be (81 ± 4)%. The 
photon identification efficiency is measured as a function 
of E T using a Monte Carlo simulation. A systematic 
uncertainty of 4% is assessed from the difference between 
the simulated electrons and electron candidates in Z — > 
ee data, and the difference between simulated electrons 
and photons. The photon identification efficiency is Ex- 
dependent and rises from about 75% at 8 GeV to about 
90% above 27 GeV. 

Backgrounds from processes where the photon is real 
and one or both of the leptons are misidentified are found 
to be negligible. Contributions from Z{— > t + t~)j events 
with leptonic decays of the tau are less than 1% of the 
sample. The only significant source of background to 
Zj production is from Z+jets processes in which a jet 
is misidentified as a photon. We estimate the Z+jets 
background by folding the jet-Ex spectrum in Z+jets 
events with the probability for a jet to be misidentified 
as a photon. The probability is measured as a function 
of the photon candidate's Et using a sample of events 
dominated by QCD multijet processes. We correct the 
misidentification probability for direct photon produc- 
tion (7+jets) by fitting the photon candidate Et distri- 
bution to the functional form derived in ^(J. For low 
Et {Et < 75 GeV) this contribution is measured to be 
9%, and we take this number as a systematic uncertainty. 
The misidentification probability is about 5 xlO -3 and 
decreases with Et- 

We use an event generator employing leading order 
(LO) QCD calculations with a detector simulation tuned 
with Z boson candidate events to calculate the accep- 
tances for the data and expected rates from both the 
SM and anomalous Zj productions ||. We use the 
CTEQ6L [HI parton distribution function (PDF) set. 
We estimate the uncertainty due to PDF choice using the 
prescription in jl2| to be 3.3%. Using a NLO Z7 Monte 
Carlo [l^l generator, we calculate an independent K- 
factor to parameterize the effect of independent NLO 
corrections in the LO Monte Carlo sample. The un- 
certainty due to the choice of if -factor (flat 
dependent) is found to be negligible. 

We observe 138 events in the electron channel, to be 
compared to the SM estimate of 95.3 ± 4.9 e + e~~j events 
and 23.6 ± 2.3 background events. In the muon channel, 
we observe 152 events vs. an estimated 126.0 ± 7.8 SM 
/i + /i~7 events and 22.4 ± 3.0 background events. The 
uncertainty in the SM signal is dominated by the uncer- 
tainty in the lepton and photon reconstruction efficien- 
cies, and that in the background estimation is dominated 
by the uncertainty in the jet misidentification probability. 

The Et spectrum for photon candidates is shown in 
Fig. n with the estimation of the total SM prediction 
and its QCD background component overlaid. The high- 
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FIG. 1: Photon candidate Et spectrum for ll'y data (solid 
circles), QCD multijet background (shaded histogram), and 
the standard model plus background (histogram). The shaded 
band is the systematic uncertainty on the SM plus back- 
ground. The Monte Carlo distribution is normalized to the 
luminosity. 



est transverse energy photon in the electron channel is 
105 GeV, while the highest transverse energy photon in 
the muon channel is 166 GeV. In Fig.[21we plot the three- 
body mass (M« 7 ) against the dilepton mass {Mu) for 
each event in the data. The ISR events with a dilep- 
ton system produced by an on-shcll Z boson populate 
a vertical band at Mu around Z boson mass, Mz, and 
Mu 7 > Mz- The on-shell Z boson FSR events clus- 
ter along a horizontal band at M« 7 = Mz and have 
Mu < Mz- The Drell-Yan events populate the diagonal 
band with Mu ~ Muj extending from the lower left to 
the upper right corner of the plot. 

For events satisfying AlZe 7 > 0.7, E T > 8 GeV, and 
Mu > 30 GeV/c 2 , the combined cross section is mea- 
sured to be 4.2 ± 0.4 (stat+sys)±0.3 (lum) pb, where 
the first uncertainty includes contributions from statis- 
tics and all systematic effects except the luminosity, and 
the second is due to the luminosity measurement uncer- 
tainty [14| . This value is in agreement with the expected 
value of 3.9lg'2 pb from NLO theory calculations [l3| . 

Given the separation exhibited in Fig. |2 we can mea- 
sure a cross section of ISR-enhanced Z7 production. 
By minimizing the effect of final state radiation, any 
anomalous contributions from a trilinear boson vertex 
would become more apparent. By requiring the dilep- 
ton mass and the three-body mass to exceed 65 GeV/c 2 
and 100 GeV/c 2 , respectively, (along with all previ- 
ous requirements), the SM Monte Carlo simulation in- 
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FIG. 2: Dilepton+photon vs. dilepton mass of Z-y candi- 
date events. Candidates in the electron channel are shown 
as empty circles, while the muon mode candidates are shown 
as stars. 



dicates that 80% of the remaining events are due to ini- 
tial state radiation. For this restricted sample we ob- 
serve 55 and 62 events in the electron and muon chan- 
nels, respectively. The cross section is measured to be 



1.07 ± 0.15(stat+sys) ± O.OJflum) pb, in agreement with 
the expected 0.94i£o5 pb 0. 

Given the good agreement observed between the data 
and the SM prediction, we extract limits on anomalous 
couplings. We generate Monte Carlo events in a two- 
dimensional grid of CP-violating anomalous couplings 
(hY and h% ) and do the same for CP-conserving (h^ 
and h\ ) anomalous couplings. We calculate the like- 
lihood of the agreement between the E^, distribution in 
data to the estimated background and Monte Carlo simu- 
lation for each point of the grid. Assuming Poisson statis- 
tics for the data and Gaussian systematic uncertainties, 
we extract the 95% C.L. limits on each of the anomalous 
couplings while assuming the others are zero. The limits 
on CP- violating and CP-conserving anomalous couplings 
are nearly identical. We also find the limits on real and 
imaginary parts of the couplings to be similar as well. 
We present the limits on both real and imaginary parts 
of the CP-conserving and CP-violating couplings in Ta- 
ble [J The two-dimensional limit contours on individual 
CP-conserving couplings are shown in Fig. [3J 

In conclusion, we have studied a sample of 290 U-y 
events, consistent with Zj production. This sam- 
ple exceeds that previously collected by D0 by an 
order of magnitude. This is due to three times 



TABLE I: Summary of the 95% C.L. limits on the anoma- 
lous couplings. Limits are set by allowing only the real or 
imaginary part of one coupling to vary; all others are fixed 
to their standard model values. As indicated, we find limits 
on CP-conserving and CP-violating parameters to be nearly 
identical. We also find that nearly identical limits apply to 
the real or imaginary parts of all couplings. 



Coupling 



A = 750 GeV 



A = 1 TeV 



\®e{ht ,3o)\, |3m(/if ,3o)l 
|Ke(/if 0i 4 )|, |Sm(ft2o,4o)| 



\3te(H 

1^(^20,40 



1,30/ 



)l, \<atm(fq 



0.24 
0.027 

0.29 
0.030 



0.23 
0.020 
0.23 
0.019 



more integrated luminosity, an increased production 
cross section associated with the 10% higher center- 
of-mass energy, and significant improvements in par- 
ticle detection efficiency achieved with the D0 Run 
II upgrade. The U-f cross section is measured to 
be 4.2 ± 0.4(stat+syst) ± 0.3(lum) pb. After addi- 
tional selection requirements, most of the final state 
radiation is removed, leaving the sample dominated 
by initial state radiation. The cross section for 
this ISR-enhanced Z7 production us measured to be 
1.07 ± 0.15(stat+syst) ± 0.07(lum) pb. These values 
are consistent with the SM expectations. We observe 
no significant deviation from the SM expectation in the 
total cross section or photon Et distribution, and thus 
extract limits on anomalous Z7 couplings. The one di- 
mensional limits at 95% C.L. for both CP-conserving and 
CP-violating couplings (both real and imaginary parts) 
are |/if 0)30 | < 0.23, |/if M0 | < 0.020, \h\ 30 \ < 0.23, and 
1^20 40 1 < 0.019 for A = 1 TeV. These limits are substan- 
tially more restrictive than previous results which have 
been presented using this formalism Q. The limits on 
hXo and h\ are more than twice as restrictive as the 
combined results of the four LEP experiments (l5| . 
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FIG. 3: The 95% C.L. two-dimensional exclusion limits for 
CP-conserving ZZ^f (a) and Z77 (b) couplings for A = 1 TeV. 
Dashed lines illustrate the unitarity constraints. 



